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Abstract 

The structure of leading nonperturbative corrections to the static Coulomb potential in 
QCD at small distances is analyzed. We argue in favor of the correction linearly dependent 
on distance and remark that lattice measurements of static potential for charges in higher 
representations can distinguish between different phenomenological models explaining its 
existence. Related problems of validity of Dirac quantization condition for running charges 
in abelian theory and significance of the quantity (A 2 ) are briefly discussed. 



1 Introduction 



Much attention has been given recently to the question about next-to-leading terms for the 
static potential in confining theory such as QCD at short distances, see [U 121 EE El OH E| and 
references therein. This problem is closely related to the structure of OPE in confining theories, 
despite to define the potential one has to go to the large time limit and leave therefore the region 
of applicability of the standard OPE 7j . In particular, the appearance of new terms in the so 
called SVZ condensate expansion jHj is discussed j2J. Phenomenologically, the static potential 
between charge and anticharge is given at small distances as the following expansion 

V(r) = — + c + c x r + c 2 r 2 + C(r 3 ) (1) 

where the coefficients c, may also depend on r but only logarithmically. Some typical patterns 
look as follows 

• Abelian charges (e; — e) in infinite space beyond tree level 

c_! = -4^ ! c 1 = c 2 = ... = (2) 

• Abelian charges (e; — e) in the cavity of the size L, tree level [5| 

e 2 

c_i = ; ci = ; a = 7;e 2 L~ ( * +1) for i > 2 (3) 

Aix 

• Nonabelian charges (gTg; —gT^) in the representation D of SU(N) interacting with soft 
nonperturbative fields JUJ 

d = -&. = -C D (g 2 /47r) ; c x = ; c 2 = ^a s (G%G%)T g (4) 



where in the last case {G a ^ u G a ^ u ) is nonperturbative gluon condensate [H] and nonperturbative 
correlation length T g characterizes the fall-off of the gauge-invariant two-point correlator of the 
field strength tensors normalized to the gluon condensate at the origin (see review ^3] and 
references therein). The eigenvalue of quadratic Casimir operator reads as 

C D = Ti D T a D T a D ■ Tr D l D = l (5) 

In the expressions above 7's stay for some dimensionless numerical factors whose values are 
of no relevance for our discussion. What attracts the attention in all mentioned cases is the 
absence of linear correction proportional to C\. This is also true for another models ^2] which 
we have not included in the list. 

However, one cannot propose strong theoretical arguments why c\ should vanish in the 
context of QCD or another theory with complicated vacuum structure. Moreover, there are 
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at least two independent sets of lattice data indicating that C\ 7^ [UJ [T3] . Different models 
explaining small distance linear correction have been proposed (see [1-6] and references therein). 
They can be divided into three groups. In the first group there are models (21 El E] which are 
based on different corrections to the perturbative two-point correlation function - either in the 
form of gluon propagator (A^(x) A b u (y)) or in the form of gauge- invariant path dependent field 
strength correlator (Tr G fiv (x)^(x, y)G pa (y)<&(y, x)). The second group is represented by the 
"dynamical cavity" model pQ, which essentially explores Q but with L = L(r) instead of a 
constant. Finally, studies of the static potential at small distances in abelian confining theories 
[Hj form the third group. We briefly comment on the latter one at the end of the paper, while 
our attention is focused on the first and the second groups. The main point of our remark 
is that despite many common points, it is possible to make distinction on the lattice between 
models of the first and of the second group by studying the small distance static potential for 
charges in higher representations of the gauge group. 

The paper consists of two relatively independent parts. In the first part (sections 2 and 3) 
we discuss some elements of short-distance physics in abelian theories. In particular, we present 
a short review devoted to the following question: does Dirac quantization condition for bare 
charges eg = 2tt stay intact when loop effects are included? There have been long discussion of 
this question in the literature, see ^1 and references therein. The affirmative answer is given 
and physics behind it is explained. We will also point the reader's attention in the section 3 
that local average of the square of the photon field (A ll (x)A IM (x)) can be defined in covariant 
gauge in the context of theory with boundaries. 

The second part of the paper (sections 4 and 5) is concentrated on the nonabelian theory. 
We demonstrate the difference between "dynamical cavity" — like model incorporating Meiss- 
ner effect, and "modified propagator" - like models in the section 4, which makes it possible to 
distinguish these scenarios. The section 5 presents our conclusions. 

In the rest of the paper unless explicitly stated otherwise, we work in Euclidean space with 
the notation for the four-vectors k = (hi, k 2 , k%, k±) and scalar product kp = k^p^d^. The three 
vectors are denoted as k = (ki, k 2 , k^) and the Wick rotation corresponds to the replacement 
&4 — > iko. 

2 Elements of short-distance physics in Abelian theory 
with monopoles 

Since we are interested in the physics of interactions at small distances, as a warm-up example, 
we consider here the theory of abelian vector field A^ interacting with the massive electrically 
charged matter field and external monopole currents. An integral part of this section is of 
mini-review type since most of the discussed results can be found in the literature. We believe 
however that consistent presentation of them can be of some use for the reader. 
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After integrating matter fields out, the partition function of effective low-energy theory is 
given by 

Z\fi\ = J VA.exp (-1 / d d P A,{p)[D-X{p)A»{-p) - S eff [A] + i J d d p jfa)A"(-pj) (6) 

where S e //[A] contains interaction terms for the field A^. The function [.D _1 ]£(p) can be 
decomposed as 

[D-X(P) = Kdo{p 2 ) + P^Pud^p 2 ) (7) 

where the function di(p) is gauge-dependent if the theory in question is a gauge theory. In the 
latter case the theory can be rewritten in the field-strength formulation ^ZJ as 

Z% v \ = J VF, U 5{d v F» v ) exp J d d pF^{p)^ v {p)F pa {-p)) ■ 

■ exp ^—S e ff[F] + i J d d PJ ; u (p)F^(-p^ (8) 
where F^ u = fe^/jF^ and 

A£(p) = - 6X)A (p 2 ) + (p,p% - p^K + p v p°^ - p u f5;)^{p 2 ) (9) 
The relation 

d (p 2 )=p 2 (A (p 2 )+p 2 A 1 (p 2 )) (10) 
provides correlator's matching between the two formulations. 

The field strength formulation (JHJ) allows one to introduce monopoles into the theory per- 
forming the shift 

SiPF^^SiPFr-ff) (11) 

If one may neglects S e ff[F] (see discussion below), straightforward integration gives the result- 
ing monopole partition function: 

Z\3™,3% = 0] = / VB^exp (-1 / d d pB,(p) 6 -^^ B^-p) + i J d d p j™(p)B»(-p)j 

P (12) 

For the free theory do(p 2 ) = p 2 , A (p 2 ) = 1 , Ai(p 2 ) = 0, which means that © coincides 
with ()12j) up to the interchange of electric and magnetic currents (charges), in other words pure 
photodynamics is exactly self-dual. If we switch on the interactions, the inverse propagator 
gets renormalized and runs with p 2 as do(p 2 ) = p 2 (l + II(p 2 )) where II(p 2 ) is the corresponding 
polarization operator. Let us assume for the time being that the dynamics keeps the gauge 
invariance of the theory intact, which means, in particular that 11(0) = 0. When computing 
the static potential in d = 3 + 1 dimensions, one gets formally in the electric case 

2 f d 3 p 1 

V«{r) = -e J ^ ^ exp(.pr) (13) 
4 



while corresponding expression in the magnetic case reads 

d 3 p d (p 2 ) 



V m rn(r) = -g / — — 77 t exp(ipr) (14) 



(2tt) 3 p 4 

where the charges e,g are physical renormalized (but not running) charges and r = |r|. 

It is seen from ()13|) . (fLf)) that the leading corrections to the Coulomb potential have 
different signs in electric and magnetic cases: 

\ 5V eS (r) = -1 6V mffl (r) = ( ^± exp(zpr) (15) 
e z g A J {zttj 6 p z 

Sometimes this is considered as a proof 1 of the fact that Dirac quantization condition 
eg = 2irn holds also for running charges (see discussion and further references in ^H]) if it does 
for the bare charges: 

e go = eg = e(r)g(r) = 2vm (16) 

where eo, e, e(r) denote the bare, the renormalized and the running charge, respectively. The 
situation is much more subtle, however. To be more concrete let us take the best studied 
case of spinor electrodynamics as an example. The problem is that introducing monopoles by 
fTTf one sets the scale of magnetic fields by the large coupling g ~ 1/e, since j™ contains 
factor g ^> 1. As an immediate result, all higher order irreducible correlators of magnetic fields 
interacting with the electrically charged matter field described by S e ff can play a role. One 
should distinguish three different kinematical regions. For large distances rm ^> 1, where m 
is electron mass the corrections to the static potential are exponentially damped with r for 
both V eg (r) and V mi n(r). Of main interest is the small distance region rm <C 1. The leading 
term in the fc-point irreducible vertex in S e ff contributing to the monopole potential (fLfj) at 
small distances contains the factor g k e k = (2nn) k while the analogous contribution to (fTHj) is 
suppressed as e k e k = e 2k 1. On the other hand the k > 4 correlators describing the processes 
of multi-photon scattering do not have 2 contributions logarithmically rising with p 2 . It means 
that for large p 2 and not too large value of the product eg one can indeed take the one-loop 
result (fTB^) as a correct first approximation. With the product eg getting large, however, higher 
order terms in the non-Gaussian part S e ff are becoming more and more important and their 
effect might overcome the leading kinematical one- loop logarithm and spoil (fTtjj) . However, this 
is not what happens. Since the effective Lagrangian is formed in this case at the distances 
I ~ f/y/eg which are for eg ^> 1 smaller that the typical scale of the field change r (see, e.g. 
JB1EI]), the constant field approximation can be used. The Lagrangian in constant magnetic 
field is given (we use Minkowskii metric conventions here) at the first order by the following 
expression [20] 

HV l-flog^ (IT) 



37T m 2 , 

1 Actually only at the next-to-leading order in electric coupling e 2 , which is assumed to be a small parameter. 
Notice that II (p 2 ) is 0(e 2 ). 

2 As is well known, at k — 4 the amplitude formally diverges logarithmically, but this divergence is exactly 
cancelled in the sum of all diagrams 
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where 4ira = e 2 and H = |H|. Since the point-like magnetic charge is defined microscopically 
via divergence of H, div H(r) = g ■ S(r), one gets from ()17|) for the running magnetic charge 
with logarithmic accuracy 



On the other hand, point-like electric charge is defined via Maxwell's equation div D(r) = 
e • S(r) where D = dL/dE, which leads to the standard Uehling-Serber result plj (again with 
logarithmic accuracy) 



This is precisely the same answer one can get from (JT5j) with the standard one-loop expression 



At first glance the coincidence between (|15|) and (|18|h (|19p seems rather surprising since 
the latter two expression were obtained from the exact strong field Lagrangian (|17|) which sums 
up infinite number of graphs while the former one is a trivial outcome of Gaussian integrations 
with one-loop two-point function. The reason for this result is the remarkable correspondence 
between QED at small distances and QED in strong fields, studied in the series of papers 
[T9*] . It follows directly from the minimal coupling principle for the gauge and charged matter 
fields: Pn — > = p^ — eA^ which gives a hint that large p^ and large eA^ can describe 
one and the same physics. The above example shows rather nontrivial manifestation of this 
correspondence. It is worth saying that the transition from one-loop dominated regime (jl5j) to 
truly nonperturbative one (|T7|) at fixed distance r with the increasing of n — ^ is a physical 
phenomenon, not related to any formal redefinitions of coupling constants etc. Notice also that 
no topological objects like Dirac strings have been involved in our analysis. Unfortunately, up 
to the author's knowledge there still exists no rigorous proof of in QED in all orders of the 
electric coupling and/or beyond logarithmic accuracy. 

If H(p 2 ) = p 2 U(p 2 ) does not vanish when p 2 — > 0, the theory leaves Coulomb phase and 
one cannot introduce external monopoles simply by A well known example of such theory 
is given by abelian Higgs model where the electrically charged field is condensed: 



Attempts to write down exact expression for the confining potential in this theory encounter a 
serious problem. The reason for that is the physical one: since the confining string is created 
between the particles (external monopole-antimonopole pair), its quantum dynamics should be 
properly taken into account. The confining interaction can not therefore be described in terms 
of particle exchanges. Moreover, even in effective abelian theory framework there exists no 
consistent procedure to perform corresponding stringy calculations analytically at the moment. 
The best one can do is to compute the Wilson loop for a particular choice of the confining 
string world-sheet geometry. The conventional choice is the minimal surface S for the given 
contour C, the flat one if the contour is rectangular (as it is for the static potential). In the 






(19) 



for n( P 2 ). 




(20) 
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London limit the Higgs boson is much heavier that the vector boson, m# ^> m = erj and such 
defined potential reads in this limit (see, e.g. [22]): 



KlSM = 9 2 Um I / A / (Py (i E^(x)A m (x - j,)^,)) - j J^™^ (21) 
where the surface current £ M/ ,(x) and the kernel A m (x — y) are given by the expressions 

= / da^) 5 {4 \x(£) -x) ; A m (x) = / — — - — exp(ipx) (22) 

J J (27TJ 4 ;r + m z 

S 

We have indicated by superscript ^ that the potential depends on the profile of the string chosen 
in (|22j) . At large distances rm ^> 1 the potential (J21j) has linear asymptotics and describes 
confinement. 

We are to address the question of applicability of (}2*Tj) at small distances. Obviously, 
the expression ()21|) breaks down at r < m^ 1 since at such small distances the scalar field 
condensate gets excited and virtual loops of scalar particles contribute to the polarization 
of the vacuum. If one takes (|2T)|) as some low energy effective theory, the ultraviolet cut-off 
implicitly enters through higher- dimensional operators. Moreover, if the Higgs boson in this 
effective theory is composite in terms of underlying microscopic theory degrees of freedom, the 
vector bosons with the energy higher than the corresponding binding energy start to resolve the 
constituents. In any of these cases the potential (|2*T|) gets modified. Let us consider the simplest 
possible modification which is to replaces m 2 in (J21j) . (J22)) into vector boson self-energy II(p 2 ) 
such that 11(0) = m 2 . Such replacements is in line with (I13j) . (J 14)) and can be justified 3 in terms 
of perturbation theory in the electric coupling e. Needless to say that the high-momentum 
asymptotics of n(p 2 ) will be different for different microscopic scenarios outlined above. 

The leading short- distance contribution to (|2*Tj) is attractive Coulomb potential while there 
are two kinds of sub-leading corrections: 

V [ ^ T \r) = + 6v£l(r) + 5V m 2 l{r) + 0{g 2 e') (23) 

The second correction term comes from the second Yukawa-like term in (|21)) and formally reads 
as 

J^vSUr) = / exp(!pr) (24 ) 

The low - p 2 divergence is artificial and can be easily regularized. The first term in (|21|) 
produces the following correction 

L 5v m (T) = r .lfh.^Sl 

g2 ov mm (r) J (2jr)s k2 

3 Of course, 11(0) = in the perturbative QED context 



'4 sin 2 (Jf J 



(h) 2 



(25) 
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where for convenience we put the Wilson contour in the (3,4) plane. It is seen that the 
dependence of the confining term - induced correction ()25|) on r is completely determined by 
the behavior of II (p 2 ) as a function of p 2 . In particular, if II (p 2 ) is such that the r - dependence 
of the integral is weak enough, the correction to the potential will be approximately linear. If 
fl(p 2 ) = m 2 , the Yukawa-like correction is linear in r and negative: SVmfLi 1 ") = —(g 2 m 2 r)/(8ir) 
(we have omitted irrelevant constant term). This is the only correction (up to the change g — > e) 
one would get for the static electric charges without taking into account the string dynamics. 
The "confining" correction (J25|) is positive in this case, but logarithmically divergent in the 
ultraviolet region. This divergence is to be cured if the correct form of ITQo 2 ) in high-momentum 
regime is taken. Presumably, the relative minus sign between 5Vmm{ r ) and SVmm{ r ) persists in 
this case as well. 

It is worth to put an emphasis on the fact that the physics behind 5Vm-k{ r ) and 5Vm-k{ r ) 
is very different. In a sense, the former one describes the process of string tension formation, 
which starts essentially at some ultraviolet scale in this model and goes all the way up to the 
distance scale given by m~ l . Since the corresponding formation "speed" is only logarithmic, one 
can have indeed approximately linear potential at small distances - "ultraviolet confinement", 
the phenomenon which apparently goes beyond the conventional OPE. 



3 Comment on the quantity (A 2 ) 

The linear correction to the potential imply the existence of local d = 2 parameter in the 
theory proportional to C\. It cannot be local polynomial function of original fields of the 
theory since the only dimension two candidate, average of vector potential squared (Tr A^(x) 2 ) 
is not gauge-invariant. Moreover, there is a gauge, defined by infinitesimal transformation 
U{x) = 1 + igA^{xo){x — x$Y + .., such that u A®(xo) = 0, i.e. one always has locally 

man ( u Afaj) 2 } = (26) 

where G is the gauge group, for abelian or nonabelian theory. In abelian case, however, the 
minimum of integral of (^(a;) 2 ) over the whole space can be rewritten in terms of gauge- 
invariant quantities |2Sj- The physical relevance of such nonlocal object for nonabelian case is 
under debate jMl23]. 

We would like to point the reader's attention to the fact, that despite it seems to be 
impossible to assign any gauge-invariant meaning to the local quantity (A 2 ), one can have in 
some cases ^-independent definition of (A 2 ), where £ is the covariant gauge fixing parameter. It 
happens when one studies the theories which contain some external parameter of the dimension 
of mass/length and if the part of the effective action depending on this parameter does not 
depend on £. By way of example let us study photodynamics with the Casimir-type static 
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boundary conditions (see review |2Hj and references therein): 

Z[# = / VA,6[BC] exp (- J d d x (±F, v {x)F^{x) + 1 (c^(x)) 2 + (27) 

where the boundary conditions read as 

= S(n^F flv (x 1 ,x 2 ,x 3 = ai,x 4 ))5(n fJ ,F fiv (xi,x 2 ,x 3 = a 2 ,x 4 )) 

Here n M = (0,0, 1,0) is a unit vector in the ^-direction and a\ , a 2 mark the X3 - positions 
of parallel infinitely thin ideally conducting plates. The distance between plates is given by 
a — \a 2 — Oi|. 

The gauge field propagator can be readily obtained from (|2~7|). it consists of two parts: 

(^(x) = D^(ar, j/) = £>2, } (x - y; £) + ^VO, y; a u a 2 ) 

where the term Dffj(x — y;£) is the standard gauge-dependent tree-level photon propagator, 
while the term D^ix, y; a\, a 2 ) encodes the information about the boundaries. The exact form 
of the latter was first obtained in [2Zj. The function D^J(x — y; £) is translation-invariant but 
gauge-dependent, while the function D^x, y; Oi, 02) depends on x% and y% separately, but is 
^-independent (not to be confused with the gauge-invariance!). Moreover, it has finite limit 
when x approaches y. Therefore one can address the issue of the (A^ 1 (x)A fl (x)) condensate in 
the Casimir vacuum exactly in the same way one computes the energy density in this problem, 
namely, subtracting the boundary- independent part: 

(A 2 (x)) de d lim ((A^x)A»(y)) - (A„(x)A^y))^) = Dftx, x; a u a 2 ) (28) 

It is convenient to introduce the notation z = ai + a 2 — 2x3. When, using exact expression for 
D IJJ/ (x,y;a 1 ,a 2 ) one gets: 



12a 2 



3 2 ^ 7TZ 



1 + - tan z — 



(29) 



2 V2a, 

Such defined local dimension two "condensate" is ^-independent, strictly positive and diverges 
on the boundaries. In completely analogous way one can define (A 2 )-" condensate" at finite 
temperature. 

In the context of Casimir problem one is usually interested to study the changes of physical 
quantities like components of energy-momentum tensor with respect to their vacuum values. 
The above example demonstrates that "non-physical" quantity like (A 2 ) is also nontrivially 
modified. It is worth noting that (|29|) is not related to any linear correction to the potential - 
the potential for dipole between the mirrors is of the form (JHJ) with L = a. 



4 "Dynamical cavity" vs "modified propagator" 

One of the main assumptions of the SVZ framework can be formulated as a wordplay: the 
nonperturbative vacuum of the theory is not perturbed by the external sources. In case of large 
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virtualities it is sometimes justified using the language of condensed matter physics: since the 
time scale of the hard process we wish to study is much smaller than typical relaxation times 
characterizing such a medium as nonperturbative QCD vacuum, it is reasonable to suggest that 
the vacuum state remains unchanged. 

It is of interest to relax this assumption and to study the corresponding physics. An 
example of the situation where one has to do it is provided by the monopole-antimonopole 
pair in the superconductor. It is well known that the condensate must be broken along some 
line connecting the particles, whatever small their charge is (see review jTHj and references 
therein). The reason for that, eventually, is the nonperturbative nature of the interaction 
between magnetic and electric particles. 

First, we re-derive the results of "dynamical cavity" model pQ in more quantitative way. 
After that, we will come to the "modified propagator" model Since we consider gluodynam- 
ics in this section, we switch to dual terminology, so that particles which are confined (quarks) 
carry electric and not magnetic charge. The microscopic description of confinement as dual 
Meissner effect via monopole condensation refers to the abelian projection procedure [2H1- We 
adopt the physical picture of confinement as dual Meissner effect but will make no use of the 
abelian Higgs model-motivated language in this section. As it is well known the energy density 
of nonperturbative fields is given by the energy-momentum tensor trace anomaly 

-in \ P( as 

4 W 16a, 



f = T(W = ?(<v(W)> ( 30 ) 



„/ \ nda s (ii 2 ) /ll „ r 2 ^\ a 2 „. ,. . , 



where 

= n 2 - 

dfi 2 V 3 3 J J 2tt 

In principle the average (J5U|) depends on the actual state of the system, for example, one 
can study its density or temperature dependence. In terms of invariant functions D(z 2 ) and 
Di(z 2 ) the vacuum gluon condensate is given as [TU] 

(g 2r Ti GV(0)G^(0)> = 24(L>(0) + DM) (32) 

where only the function D(z 2 ) is responsible for confinement. In particular, the deconfinement 
transition is characterized by the condition D(z 2 ) = 0, with the smooth behavior of D\(z) over 
the phase transition [30 . 

One can rise a question what happens to (}30|) in color field of external source. Owing to 
SU(3) and 0(4) invariance, such average can be defined as 

r M / as r* Mr>* m\ 9 * ( Tr (gg^OOggjg)) ; Tr p ex P j Ic A P dz P ) ( oo\ 
G2{r) ~ \7 G ^ (r)G > (r) / " 2^ (TrPexp^ c ^ p ) (33) 

where C stays for the (closed) world-line of the source. If the point r is far enough from C, the 
average reaches the vacuum value: (^(r) — > G2. Correspondingly, one defines e(r) = e(E(r)) 
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via (jHOJ). To the best of author's knowledge, no systematic lattice analysis of behavior of e(r) 
in the strong external field has been performed. This problem is of direct relevance for our 
discussion. Namely, let us put a single static charge in the nonperturbative QCD vacuum. Due 
to confinement property there should be an anti-charge somewhere in the space to end the 
confining string and therefore the field distribution will not be spherically symmetric. Suppose, 
however, that we look at the gluon fields in a thin spherical layer between r and r + Ar, 
where r is smaller that the typical nonperturbative nonlocality scale given by T g (it is worth 
mentioning that the same quantity T g defines the width of the confining string). Then it is 
natural to assume that the Faraday flux lines of the charge are affected only weakly by the 
opposite charge sitting on the other end of the string. On the other hand, they are affected by 
the nonperturbative vacuum and act back on it. Correspondingly, we take the energy of the 
gluon fields as a sum of perturbative and nonperturbative parts 

£ = £ p + £ np = ~ J dh E a (r)E a (r) + J d 3 r e(r) (34) 
v v 

where E a (r) is the perturbative electric field of the charge. In principle, the coordinate de- 
pendence of E a (r) and e(r) = e(E(r)) on r can be very complex, resembling the formation of 
intermediate state in the vicinity of the charge, like it happens in ordinary superconductors. 
It can be energetically advantageous to squeeze slightly perturbative flux lines at one part of 
the volume 4 and to rarefy them in the rest of it (this increases the first term in (}34|) ) but com- 
pensate the excess of energy by the nonperturbative energy gain from the rest of the volume, 
which can be possible if e(r) decreases towards its vacuum value when E a (r) decreases in this 
part of the volume. This is the essence of Meissner effect. Let us simplify the matter further 
by approximating as 



-4 



C 



D 
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4yrr 2 \ 2 



Mr), 



■ S^Ar + ex ■ 3 1 (r)Ar + e 2 • S 2 (r)Ar + 0(Ar 2 ) (35) 



We have taken e(r) as constants ei and e 2 in the parts of the volume SiAr and .S^Ar, respec- 
tively. Notice that S\(r) + S%(r) = Anr 2 . The expression fl33j) assumes that the parts of the 
space occupied by the nonperturbative fields with the (vacuum) energy density e 2 completely 
expel the perturbative field flux lines into the regions where the energy density of nonperturba- 
tive fields is given by ej. The rescaling of electric perturbative field of the point charge reflect 
the flux conservation. 

Since S\(r) can be never greater than 47rr 2 , there exists critical radius (corresponding to 
the stationary point of (jUSI) 

< = . *' ; (36) 
87r(ei — e 2 ) 

When r reaches r c from above, the region of unperturbed vacuum S 2 shrinks to zero. Physically 
the model describes the formation of a bag of the radius r c , there perturbative fields are strong 

4 Since we are working at small distances, the total perturbative flux (i.e. the number of the flux lines) is 
conserved. More accurately, because of asymptotic freedom it is not, but its violation due to asymptotic freedom 
is weak logarithmic effect 



11 



enough to change the nonperturbative vacuum state significantly. We are not addressing an 
interesting question at the moment, is there local deconfinement transition (in any sense of 
the word) in this volume (what is natural to think of being based on abelian Meissner effect 
analogy). Let us estimate r c numerically. We use different sets of data [3T] and take for the 
gluon condensate 

G 2 = /^<™,G« \ = (0.014 4- 0.026)GeV 4 (37) 

Hence for the nonperturbative vacuum energy in one loop choosing N = 3, Nf = 2 we obtain 

ei - e 2 = «e - e = (1 - «) • (4 4 8) ■ 10~ 3 GeV 4 (38) 

The value of the positive factor « = e\je 2 < 1 is unknown and enters as free parameter. Taking 
a s = 0.47 (corresponding to the charges in fundamental representation and a s (M T )), we finally 
get 

r c = (1 - . (0.2 0.3) Fm (39) 

Unless k is unnaturally close to unity this rather small value of r c is compatible with the 
discussed picture - had it been much larger than T g , the analysis would be meaningless. 

Let us, following [1 , consider now a small color dipole of the size R. The expression for 
the energy (|34|) stays intact, but the physics of critical distance r c is different. Since there 
is no flux of perturbative field through the surface of the volume with the dipole inside, it is 
advantageous to confine flux lines inside some finite volume. Suppose that we have chosen the 
volume of some fixed size L > R around the dipole, such that vacuum energy density e 2 = e 
outside is reduced to e\ > e 2 inside it. The perturbative electric field felt on the boundary of 
this volume is proportional to R/L 3 . It is important that it goes to zero when R decreases. 
It means that making R small enough the total energy can also be lowered by appropriate 
decreasing of L, i.e. L = L(R). Indeed, when we go from L to L — AL the perturbative energy 
increase due to the rearrangement of perturbative field is given by 

A£ * = \ c ° (i) 2 1 d ' r (* + 3 ^Pl = & 'ij AL + ° (AL2) (40) 

L-AL ^ ' 

while the nonperturbative energy gain is 

AS np = 4nL 2 (e 1 -e 2 )AL + 0(AL 2 ) (41) 
This defines the stationary point 



aM 2 V /6 



V47r(ei-e 2 )y 

We can now estimate the nonperturbative correction to the energy of the dipole as 

R 2 , , 4vr , 4 



(42) 



8£{R) « a a ^3 + (ei - e 2 ) • ^-L 3 C = j ■ yjnafa - e 2 ) • R (43) 
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This correction is linear in R, as one can alternatively see inserting ()42j) into (f3*J). Numerically 

we get 

- • ^7ra s (ei - e 2 ) = (1 - k)s • (0.10 - 0.14) GeV 2 (44) 

The discussed picture reminds the old MIT bag model [22] ■ This similarity is formal, in 
some sense. While in the MIT bag model the bag has the typical hadron size, the "dynamical 
cavity" of the size L C (R) given by (J42|) is an object, which has physical meaning only at the very 
small distances. Since L c ~ i? 1//3 , at small enough R we are always in T g > L c > R regime, 
where all the picture is meaningful. With R and L c rising and reaching T g , the confining string 
formation process starts and we leave the domain of qualitative applicability of the model. 

We can now come to the "modified propagator" models [21 El E] Actually, their analysis 
is simpler and we take as an example the model of [2] . It is suggested to modify the one-gluon 
exchange propagator by adding "tachyon" mass term to the gluon, i.e. (in Minkowskii metric 
and Feynman gauge) 

(A a „(k)Ai(-k)) = -> V , u 5 ab (1 + ^ (45) 

Notice that (|45|) is large momentum expansion so one actually never gets closer to the artificial 
"tachyon pole". The analysis of existing phenomenology leads to A 2 ~ —0.5 GeV 2 estimate for 
the "tachyon mass" [2]. However the physical origin of the A 2 -term in the context of QCD has 
not been clarified. Alternative scenarios suggest fi 2 /z 2 - contribution 5 to the function D\(z 2 ) 
jS], strong coupling constant freezing jB], infrared renormalons pQ. It is easy to see that all 
these proposals lead to the Casimir scaling law for the corresponding correction, for example, 
with expression (|43|) one gets 

(46) 

/l A 2 r\ 

V D (r) + const = -a s .(- + _!+ Q(r 2 ) (47) 

where a s = C D a s = C D (g 2 /47r), i.e. V D (r) ~ C D . 

The magnitude of corrections are close to each other in both models, for example taking 
value (jUJ) one gets in terms of (f43j) A 2 = — (1 — /c)a • (0.4 -4- 0.6) GeV 2 which is surprisingly 
close to the value found in [2 . The important difference, however is different Co- dependence. 
The remarkable property of ([4H|) is its square root dependence on the eigenvalue of Casimir 
operator Co- This result is known in the context of MIT bag model jSS], where it takes place 
for the slope of confining linear potential. As such, square root law was definitely ruled out by 
precise recent studies of static potential at large distances in different representations of the 
gauge group (HUES]- Instead, the Casimir scaling phenomenon (proportionality of Vd{t) to 
Cd) was confirmed (see jHHl for review). However there exist no calculations of static potential 

5 It can be shown that l/z 2 term in the function D(z 2 ) would produce ultraviolet divergence on the world- 
sheet of the confining string, while such term in D\{z 2 ) is safe. 
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for charges in higher representations at small distances . Such lattice simulation will confront 
in nontrivial way (J43j) and (|47|) . 

From general point of view the Casimir scaling and y/Co law for the leading small distance 
correction would correspond to rather different physical pictures. Both cases physically describe 
perturbative-nonperturbative interference and require modification of the standard " condensate 
ideology" in QCD. However in the latter case this modification is to be by far more radical. 
Roughly speaking, the former case corresponds to a new kind of nonperturbative corrections 
to the perturbative propagator, which are "harder" than those studied before, while the latter 
case indicates that the nonperturbative vacuum structure itself is strongly affected by the 
perturbative field when it reaches some critical value. 

To be self-contained, let us briefly mention the analysis of classical abelian Higgs model 
on the shortest distances, performed in The small distance correction to the potential was 
found numerically and fitted by linear function with a good accuracy. This behavior was linked 
to the Dirac veto, i.e. the condition for the charged condensate <ft to vanish along some line (the 
Dirac string) between the dual charges. The problem mentioned in the section 2 still persists in 
this case, because the position of the Dirac string can be chosen in arbitrary way and the answer 
for the short distance potential depends on this choice. Moreover, it is not quite clear what the 
Dirac string of effective abelian theory does correspond to in terms of the original nonabelian 
theory. In "dynamical cavity" model the linear behavior (|4~H|) also results from the Meissner 
effect but without any special role of singularities like Dirac strings. Roughly speaking, the 
effect in [3| comes from the infinitely thin line between the charges, while in the discussed case 
all the volume of the size L r contributes to it. 



5 Conclusion 

There are a few models in the literature aimed to reproduce the linear dependence with distance 
of the leading nonperturbative correction to the static potential at small distances. Comparing 
different models we call the reader's attention to the fact, that possible way to put them on test 
is to check the dependence of this correction on the quadratic Casimir for higher representations 
of the charges. Such computation can shed new light on strong interaction physics and the 
problem of confinement. 
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